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PREFACE 

This  Report  represents  Cha  culmination  of  an  affore  bagun  aavaral 
yaara  ago.  Tha  United  Statea  Air  *orca  was  examining  the  feaalbllity 
of  altarnativa  ayatama  of  mlaaila  baaing,  and  baeama  intaraated  in  accu¬ 
rately  calculating  tha  preaaurea  and  motlona  following  a  nuclear  burst 
in  or  near  water.  Such  calculations  are  naceaaary  to  eatabllah  esti- 
matea  of  damage  to  nearby  atructuraa. 

Tha  material  in  tha  Report  la  highly  technical,  and  la  baaic  to 
an  underatanding  of  tha  phyalcal  behavior  of  water  at  high  temperaturea 
and  preaaurea.  The  reaulta  praaented  should  have  lasting  utility,  not 
only  for  their  original  purpose  of  damage  assessment,  but  in  a  wide 
variety  of  applications. 

This  Report  is  the  third  and  last  In  a  sariea  that  documents  re¬ 
sults  on  the  determination  of  thermodynamic  relations  for  water;  the 
others  are  RM-4969  and  RM-5050. 
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SUMMAR7 

The  BKW  equation  Is  an  empirical  relation  for  the  pressure-volume- 
temperature  equation  of  state  of  a  nonideal  mixture  of  gases.  This 
equation  and  some  of  its  nunerous  modifications  have  been  used  as  a 
basis  for  determining  thermodynamic  relationships  for  gaseous  mixtures 
at  high  pressures  and  densities.  In  this  Report,  we  test  the  BKW  equa¬ 
tion  and  some  of  its  recent  modifications  for  use  in  determining  a  set 
of  thermodynamic  relations  for  water  over  a  wide  range  of  pressures, 
temperatures,  and  densities.  On  the  basis  of  the  results  of  these 
tests,  we  propose  a  further  group  of  modifications  that  would  provide 
improved  thermodynamic  relations  for  use  in  calculations  simulating 
nuclear  bursts  in  or  near  water. 

We  test  our  various  proposed  modifications,  and  recommend  one  of 
them  on  the  basis  of  our  results.  Finally,  we  calculate  a  pressure- 
specific  volume-internal  energy  relation,  based  on  our  modified  form 
of  the  BKW  equation,  that  can  be  used  above  approximately  300  kllobars. 

The  detailed  results  of  the  study  are  summarized  on  pp.  50-52  of 
this  Report. 
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SYHBOLS 

A  ■  Halmholz  free  energy  of  mixture 

e  "  number  of  epeeles  in  mixture 

E  ■  intemel  energy  of  mixture 

T  "  Gibbs  free  energy  of  mixture 

H  ■  enthalpy  of  mixture 

h  "  specific  enthalpy 

M  -  mass  of  one  gram-mol  of  water 
o 

n  “  total  tools  of  mixture 

g 

n^  •  mols  of  the  i**1  component  in  mixture 
p  -  pressure 

p°  “  standard  state  pressure  (1  atm) 

pQ  «  pressure  at  center  point  of  the  Hugoniot  (1  atm) 

R  *  molar  gas  constant 

r  ■  distance  between  centers  of  two  molecules 

S  “  entropy  of  mixture 

T  ■  temperature 

T  -  298.16’K 
o 

V  ■  total  volume  -  M  v 

o 

V°  ■  total  volume  of  mixture  at  pressure  p°,  V°  ■  n^RT/p° 
v  -  specific  volume 
x^  "  tol  fraction  of  itl*  species 
y^  ■  mols  of  atomic  hydrogen 
y j  "  mole  of  atomic  oxygen 
y^  •  mols  of  electrons 
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chemlcal  potential  of  j^5  species  in  mixture 
stoichiometric  coefficient  of  j6*1  species  in  iC^  reaction 
lntermolecular  potential 


Superscripts 

'  •  differentiation;  e.g.,  <e*(T)  ■  dic/dT 

~  -  molar  quantity 

o  -  material  at  temperature  T  in  its  standard  state  (ideal 
gas  at  1  atm) 

*  ■  a  residual  (Imperfection)  contribution  to  a  thermodynamic 
quantity  (except  for  T*) 

Subscripts 

o  *  ambient  conditions  -  liquid  water  at  298.16*  and  1  atm 

■  index  representing  (i,j,k)ttl  reaction,  or  (i.j.k)6*1 
species 


J> 


I.  INTRODUCTION 


The  material  presented  here  represents  results  from  a  portion 
of  the  studies  that  have  been  directed  toward  the  development  of  equa¬ 
tions  of  state  of  water  suitable  for  high-energy  and  hydrodynamic  com¬ 
putations.  In  this  report,  we  review  current  forms  of  the  BKW  equation 
of  state  (discussed  in  Che  following  section)  and  recemmtend  some  modi¬ 
fications  of  its  currently  accepted  form.  We  also  test  both  the  exist¬ 
ing  BKW  equation  and  our  modified  versions  for  their  suitability  in 
hydrodynamic  computations  over  a  wide  range  of  pressures  and  temperatures. 
Although  these  results  were  obtained  some  time  ago,  they  are  still  rele¬ 
vant  to  the  major  unsolved  problem  of  obtaining  the  equations  of  state 
needed  to  theoretically  predict  the  motions  and  stresses  following  a 
nuclear  burst  in  water.  Results  based  on  the  studies  presented  here 
cover  a  significant  portion  of  the  entire  pressure  and  temperature 
range  needed  in  such  computations.  A  recommended  set  of  equations 
covering  the  entire  range  of  interest  and  containing  the  results  of 
this  study  will  be  published  in  the  near  future. 

The  following  investigation  is  based  on  the  thermochemlcal  pro¬ 
gram  developed  in  Ref.  1.  Using  the  general  outline  of  that  reference, 
we  divided  the  pressure-volume  plane  into  a  series  of  distinct  regions, 
shown  in  Fig.  1.  In  the  hlghest-pressure  region,  we  used  the  equation- 
of-state  results  of  A.  Latter  and  R.  Latter  based  on  the  Thomas-Fermi 
method.  Below  the  Thotnas-Fermi  region,  the  molecule  has  decomposed 
into  the  twenty-three  species  of  Table  1,  which  are  assumed  to  be  in 
chemical  equilibrium  for  temperatures  above  500Q*K  (Region  I  of  Fig.  1). 
Calculations  using  the  BKW  equation  of  state  have  shown  that  more  than 
99  mol  percent  of  the  mixture  consists  of  six  species  (H,  0,  H^,  O^, 

HjO,  and  OH)  at  temperatures  of  5000*K  and  7000*K.  Therefore,  in  Re¬ 
gion  II  of  Fig.  1,  only  these  six  species  are  assumed  to  exist.  Below 
2000*K  (Region  III) ,  BKW  calculations  show  that  the  H2O  molecules  re¬ 
main  essentially  intact.  Thus  the  2000*K  isotherm  is  established  as 
the  lower  boundary  of  Region  II.  Regions  I  and  II  cover  the  greatest 
portion  of  the  pressure-volume  plane  over  which  the  BKW  equation  of 
state  will  be  modified  and  tested  to  determine  thermodynamic  relations 


Pressure  (megabars) 
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Fig. )— Approximate  regions  of  the  BKW  equilibium  studies.  All  values 
Shown  of  the  sketch  are  approximate.  Isotherms  and  Hugoniot  shown 
are  illustrative  only. 
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SPECIES  AND  REACTIONS  FOR  THE  CHEMICAL  EQUILIBRIUM  PROBLEM  IN  WATER 
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for  water.  The  various  species  which  may  occur  In  Region  X  are  shown 

In  Table  1.  Stoichiometric  coefficients  for  reactions  among  the  species 

of  Table  1  are  shown  In  Tables  2  and  3.  The  upper  limit  on  Region  X  la 

roughly  established  at  20  megabars,  corresponding  to  our  estimate  of  the 

lower  limit  of  validity  of  the  Thomas-Feral  calculations  of  A.  Latter 
/ 2 ) 

and  R.  Latter'  for  H^O  at  extreme  pressures. 

The  reasons  for  terminating  the  low-pressure  results  of  Latter  and 

Latter  are  indicated  in  Fig.  2,  which  contains  the  Hugoniot  results  of 

Raf.  2  at  the  high  pressures  and  experimental  Hugoniot  points  at  the 
(3  4) 

lower  pressures.'  *  At  pressures  below  roughly  20  megabars,  the 

Thomas-Feral  model  shows  a  tendency  to  predict  pressures  which  are  too 

high.  This  deviation  of  the  Hugoniot  la  consistent  with  the  oft-repeated 

statement  that  the  Thomas-Feral  model  should  not  be  used  below  a  pressure 

of  approximately  10  megabars. Since  this  comment  applies  to  the 

heavier  elements,  which  more  generally  satisfy  the  assumption  of  a 

Thomas-Feral  analysis  than  the  ^0  molecule,  the  apparent  deviations 

of  Fig.  2  are  not  surprising.  (At  a  temperature  of  298*K  and  a  volume 

of  1  cm^/gm,  which  correspond  to  normal  liquid-water  conditions,  the 

Thomas-Feral  model  predicts  a  pressure  of  approximately  350  kllobars.) 

(4) 

Rice  and  Walsh  have  developed  an  equation  of  state  of  water  to  250 
kllobars  which  is  applicable  to  a  region  near  the  Hugoniot,  and  Walker 
and  Sternberg have  constructed  a  pressure-volune-energy  equation 
for  water  in  the  region  bounded  on  the  left  in  a  p-v  diagram  by  the 
Hugoniot  and  on  the  right  by  the  adlabat  which  intersects  the  Hugo¬ 
niot  at  250  kllobars.  A  p-v-T  relation  to  1  megabar  and  10,000‘C  has 

(  8) 

been  developed  by  Howard.  These  investigations  cover  only  a  small 
portion  of  the  region  over  which  we  are  seeking  equilibrium  solutions. 

To  our  knowledge,  the  investigations  of  Snay  and  coworkers ^  represent 
the  most  systematic  attempt  to  obtain  equilibrium  relations  for  water 
in  the  regions  of  Fig.  1. 

The  principal  analytical  tool  used  in  their  investigations,  the 

++ 

Becker-Kistiakowsky-Wilson  (BKW)  equation,  has  since  been  modified 
+See  Refs.  9,  10,  11,  and  12. 

^This  same  relation  is  often  referred  to  as  the  Halford-Kistia- 
kowsky-Wilson  (HKW)  equation. 


Table  2 


Note:  A  blank  apace  In  the  table  algnlfies  a  zero  value  for  the 
etolchloaetrlc  coefficient. 


and  improved  by  various  investigators  at  Los  Alamos. (13,14,15) 
attempt  here  to  capitalize  on  these  improvements  to  obtain  perhaps  more 
accurate  thermodynamic  relations  than  those  reported  for  water  in  Refs. 
11  and  12.^  We  also  attempt  to  estimate  the  range  of  validity  of  re¬ 
sults  obtained  by  our  version  of  the  BKW  equation. 

t 

In  Ref  12,  the  BKW  modifications  of  Ref.  IS  are  recommended.  As 
a  starting  point,  we  will  use  the  most  recent  BKW  parameters  of  Ref. 

14,  which  show  better  agreement  with  experimental  Hugoniot  data  than 
those  of  Ref.  15. 
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II.  THE  BECKER-KISTIAKOWSKY -WILSON  EQUATION 


For  M  grams  (n  mols)  of  a  pure  substance,  the  BKW  equation,  with  its 
o  g 

most  recent  modifications,  takes  the  form 


-  1  +  xe 


Bx 


(1) 


n  RT 

g 


ick 

V(T  +  9)‘ 


(2) 


where  the  parameters  B,  tc,  a,  8,  and  k  are  constants.  The  total  vol- 
ume,  V,  may  also  be  expressed  as  V  -  Mqv  whenever  we  wish  to  present 
results  in  terms  of  specific  volume,  v. 

The  variable  x  in  Eq.  (2)  bears  a  resemblance  to  a  familiar  ratio 
in  the  theory  of  imperfect  gases: 

Covolume^  of  M  mass  of  material 
o 


Total  volume  of  M  mass  of  material 
o 


If  the  parameter  9  were  not  present  (and  it  wasn't  in  the  original  form 
of  the  equation) ,  the  temperature-dependence  of  x  would  be  that  achiev¬ 
ed  by  molecules  with  a  repulsive  pair  potential  of  exponent  3/a: 

$(r)  *  r  a,  where  r  is  the  intemuclear  distance  between  two  inter¬ 
acting  molecules. 

For  a  mixture  of  c  species  of  composition  n  ■  {n. ,  ...,  n  }  (n  , 

~  A  C  J 

where  J  -  l...c,  represents  the  number  of  mols  of  the  j  species  in 
a  mixture  of  mass  Mq) ,  the  covolume  factor,  k  of  Eq.  (2)  is  replaced 
by  an  average  covolume  factor 


The  resolution  of  the  constant  in  the  numerator  of  Eq.  (2)  into 
a  product  will  become  clearer  when  we  discuss  the  form  of  the  equation 
for  mixtures. 

^The  covolume  of  a  gas  is  equal  to  4  times  the  volume  of  molecules 
treated  as  rigid  spheres,  and  roughly  represents  the  volume  unavailable 
to  the  centers  of  the  molecules.  The  covolume  factor  (p.  10  and  follow¬ 
ing  text)  is  an  empirical  constant  that  resembles  the  covolume. 


•10 


c 

k  "  l  ni*ti 

i-1  1  1 


where  Is  Che  covolume  factor  of  the  species. 

Equations  (1)  and  (3)  form  the  starting  point  of  our  equilibrium 
calculations: 


x 


*  I 

i-1 

v(t  +  e)° 


(3) 


The  parameters  k,  8,  a,  8,  and  all  of  the  covolume  factors,  ,  are 
regarded  as  empirical  material  constants  to  be  determined  by  experi¬ 
ments.  These  are  usually  detonation. experiments,  performed  with  a 
variety  of  explosives  with  different  product  compositions,  where  the 
above  parameters  are  chosen  to  best  predict  experimental  curves  of 
detonation  velocity  versus  loading  density,  and  Chapman- Jouguet  temp¬ 
eratures  and  pressures.^ 

The  strongly  empirical  nature  of  the  BKW  equation  has  been  the 
principal  source  of  a  number  of  warnings  Issued  against  its  use 

in  an  extended  extrapolation  of  experimental  results  either  to  high 
pressures  and  densities  or  to  materials  (such  as,  for  example,  0+^  in 
Table  1)  that  did  not  appear  in  the  experiments  used  to  choose  the 
parameters  of  the  problem.  However,  the  studies  of  Refs.  18,  19,  20, 
and  21  provide  little  motivation  to  use  more  complicated  equations  of 
state,  more  sophisticated  ir.termolecular  potentials,  or  more  physically 
plausible  mixture  theories  for  calculations  in  our  high-density,  high- 
pressure  region,  since  such  methods  have  not  led  to  results  that  were 
clearly  superior  to  those  obtained  via  the  BKW  equation  at  experiment¬ 
ally  attainable  pressures  and  densities.  Molecular  theories  for  equi- 

^For  further  information  on  this  subject,  we  refer  the  reader  to 
Refs.  13,  14,  15,  and  16.  Chapman- Jouguet  theory  is  discussed  in  Ref. 
17.  The  use  of  detonation  theory  and  detonation  experiments  to  set 
empirical  constants  in  equations  of  state  based  on  intermolecular  po¬ 
tentials  are  discussed  in  Refs.  18,  19,  20,  and  21. 
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librium  calculations  In  mixtures  below  the  Thomas -Fermi  region  ere  based 

on  the  use  of  spherical,  additive  pair  potentials  to  represent  the  forces 

acting  on  a  molecule.  Such  calculations  are  all  semicmplrlcal,  for  the 

pair  potentials  on  which  they  are  based  contain  force  constants,  which 

oust  be  established  by  experiment  for  most  materials.  In  addition,  the 

assumption  of  spherical,  additive  pair  potentials  becomes  unreliable  at 

pressures  as  high  as  those  contemplated  in  our  calculations,  even  for 

(21  22  23) 

such  normally  well-behaved  materials  as  the  noble  gases, '  *  *  and 

we  are  by  no  means  dealing  with  well-behaved  materials.  Thus,  the  phys¬ 
ical  theories  based  on  molecular  models  contain  many  of  the  undesirable 
traits  that  plague  the  BKW  equation.  It  appears  that  the  region  between 
experimental  data  and  the  Thomas-Ferml  model  cannot  be  treated  in  at.  en¬ 
tirely  satisfactory  manner,  and  that  each  of  the  current  methods  of  an¬ 
alysis  assumes  more  the  aspect  of  an  interpolation  than  a  physical  theo¬ 
ry  when  applied  to  this  region.  At  the  same  time  however,  extrapolations 
of  the  BKW  equation  to  infinite  pressures,  such  as  appear  in  Refs.  11 
and  12,  place  considerably  greater  faith  in  the  validity  of  this  equa¬ 
tion  than  it  deserves.  The  approach  taken  here  represents  a  compromise 
between  the  two  conflicting  views  and  usages  of  the  BKW  equation  at  ex¬ 
treme  pressures. 

We  tentatively  accept  the  BKW  equation  as  an  interpolation  formula 
between  1  and  20  megabars.  Its  parameters  are  then  adjusted  so  that 
its  pressure-volume  Hugonlot  curve  Joins  the  Thomas-Ferml  Hugoniot  at 
about  20  megabars,  and  the  experimental  data  at  about  one  megabar.  We 
next  test  the  resulting  form  of  the  BKW  equation  on  isotherms  and  iso¬ 
metrics  to  see  how  well  our  modifications  reproduce  known  results  at 
high  and  low  pressures.  We  then  recoomand  use  of  the  unmodified  or 
modified  forms  of  the  BKW  relations  only  for  those  thermodynamic  vari¬ 
ables  that  show  reasonabls  agreement  with  known  results,  and  only  over 
the  estimated  ranges  whers  agreement  does  occur.  As  we  will  show  later, 
the  Hugoniot-matching  conditions  require  treating  certain  BKW  param¬ 
eters,  hitherto  held  constant,  as  functions  of  one  or  more  thermody¬ 
namic  coordinates. 
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III.  VARIABLE  BKW  PARAMETERS 

To  decide  which  of  Che  parameters  k,  8,  a,  8,  and  k^  should  assume 
Che  role  of  variables  rather  Chan  constants,  and  what  their  functional 
forms  should  be,  we  use  the  following  three  guidelines:  (1)  We  will 
seek  simple  working  equations.  (2)  We  wish  to  vary  as  few  parameters 
as  possible  and  still  achieve  our  previously  stated  objective  of  forcing 
BKW  results  into  high-  and  low-pressure  results  of  other  solutions  or 
data.  The  more  terms  we  vary,  the  more  we  defeat  our  first  goal  of 
maintaining  simple  forms  for  computation.  (3)  We  prefer  to  vary  those 
parameters  which  lend  themselves  to  some  physical  interpretation.  In 
this  way  we  can  more  confidently  predict  the  general  consequences  of 
an  assumed  variation.  A  physical  interpretation  of  the  parameter  usually 
places  reasonable  bounds  on  the  range  of  values  which  may  be  assumed 
by  that  parameter.  For  instance,  we  know  that  the  covolume  factors 
should  always  remain  positive.  Therefore,  a  range  of  <  which  includes 
negative  values  would  be  unreasonable.  It  would,  on  the  other  hand, 
be  difficult  physically  to  justify  the  prohibition  of  negative  values 
for  the  range  of  the  parameters  6  and  8  for  temperatures  larger  than  8. 

The  individual  terms  *k  are  proportional  to  the  molecular  volume 

th  i 

of  the  i  component  of  the  mixture.  Each  of  these  individual  terms 
could  be  varied,  but  there  are  twenty-three  in  all,  and  assuming  an 
independent  variation  of  each  would  create  computational  difficulties 
well  beyond  those  that  we  would  care  to  encounter.  Besides,  there  is 
little  information  available  to  us  to  decide  how  each  species  should 
change  relative  to  all  of  the  others.  The  common  factor,  k,  on  the 
other  hand,  can  be  varied  without  introducing  too  many  subsequent  dif¬ 
ficulties  into  our  computer  programs,  and  can  still  be  thought  of  as 
yielding  a  physical  effect.  We  can  consider  the  molecules  of  the  mix¬ 
ture  as  becoming  "softer"  and  "harder,"  respectively,  with  decreasing 
and  increasing  k.  The  variation  of  k,  rather  than  the  individual  molec¬ 
ular  constants  k^,  implies  that  the  molecular  volume  of  each  species 
grows  or  decreases  proportionately,  which,  though  not  entirely  accurate, 
represents  at  least  a  plausible  manner  for  molecules  of  a  mixture  to 
change  their  size.  The  parameter  8  has  little  physical  significance, 
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except  what  can  be  gleaned  froa  Eq.  (1)  and  a  knowledge  of  the  fact 
that  x  2  0.  We  see  froa  this  equation  that  an  Increase  in  S  represents 

mm  mm 

an  increase  in  the  compressibility  factor,  pV/HX,  and  thus  a  greater 

M  mm 

departure  from  ideality  (i.e.,  the  ideal  gas  relation,  pV  -  KT) .  Treat¬ 
ing  8  aa  a  variable  would  not  severely  complicate  the  forms  of  the  equa¬ 
tions,  for  it  does  not  appear  in  the  definition  of  x.  The  parameter  a 
appears  in  the  BKW  equation  in  a  considerably  aor:  complicated  way  than 
8  or  x.  We  therefore  do  not  care  to  vary  a  if  we  can  avoid  it.  Fin¬ 
ally,  the  parameter  9  cannot  be  assigned  any  kind  of  physical  interpre¬ 
tation  that  we  can  think  of,  and  in  addition  would  have  to  vary  strongly 
at  the  higher  temperatures  to  have  any  significant  effect  on  the  solu¬ 
tions,  We  therefore  hold  9  constant  at  its  current  value  of  400*K. 

On  the  basis  of  these  remarks,  we  have  elected  to  vary  the  parameters 
k  and  6  and  to  treat  k^,  a,  and  9  as  constants. 

There  still  remains  the  problem  of  choosing  the  independent  vari¬ 
ables  for  the  functions  defining  k  and  B.  The  independent  variables 
of  the  entire  equilibrium  problem  are  temperature  and  volume.  These 
would  appear  to  be  the  best  choice  of  independent  variables  for  k  and  8. 
The  residual  internal  energy  and  residual  Helmholz  free  energy  are  ob¬ 
tained  by  integrating  the  expressions  T(3p/ST)  -  p  and  p,  respectively, 

v  t 

with  volume  (v)  at  constant  temperature  (T)  and  composition  n.  This 
integration  is  most  efficiently  carried  out  for  the  BKW  equation  by 
replacing  V  by  x  as  an  independent  variable  via  the  differencial  re¬ 
lation 

dx 
x 

T,n  «  constant 

'rnm 

The  form a  of  the  functions  would  therefore  be  subject  of  simpli¬ 
fication  if  <  and  8  were  functions  of  x  rather  than  V.  In  addition, 
the  value  of  x  is  a  better  index  of  the  departure  of  the  material  from 
ideality  than  V.  Indeed,  if  6  were  not  a  function  of  T  (and  it  will 
not  be  in  our  final  results) ,  x  would  be  uniquely  related  to  the 

^Sae  Ref.  1  for  a  discussion  of  these  integrations. 
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compressibility  factor,  pV/RT.  Therefore  T  and  x  appear  suitable  choices 
for  the  independent  variables. 

Although  we  have  narrowed  the  choice  down  to  the  two  functions 
K  -  k(x,T),  8  -  8(x,T) , 

these  are  still  too  general  for  the  ntmerlcal  study  needed  to  establish 
their  approximate  forms.  Once  again,  in  the  name  of  convenience,  we 
limit  ourselves  to  two  simple  possibilities: 


B  -  8(T) ; 

8  -  tc(T) 

(4a) 

8  ■  6(T) ; 

K  ■  K(x)  . 

(4b) 

The  equilibrium  relations  for  the  pressure.  Internal  energy,  Helm- 
holz  free  energy,  chemical  potentials,  and  equilibrium  constants  as 
functions  of  temperature  and  specific  volume  will  be  presented  in  the 
next  section  for  relations  (4a)  and  (4b) . 


15 


IV.  SUMMARY  OF  WORKING  EQUATIONS 

Th«  thermodynamic  equations  developed  in  this  section  pertain  to 
a  mixture  composed  of  the  twenty-three  species  appearing  in  Table  1, 
which  are  in  chemical,  thermal,  and  mechanical  equilibrium,  and  whose 
parameters  fall  into  one  of  the  two  following  categories: 

A:  a,  0,  k^  ■  constant 

8  -  8(T) ,  x  -  ie(T) ; 

B:  a,  0,  k^  ■  constant 

8  -  B(T),  k  -  x(x). 

Because  the  equations  for  Region  II  and  the  equations  correspond¬ 
ing  to  constants  8  and  k  are  all  special  cases  of  the  relations  which 
we  will  obtain,  there  is  no  need  to  diacuss  these  cases  separately. 

The  equations  in  terms  of  an  arbitrary  p-v-T  relation  are  devel¬ 
oped  in  Ref.  1.  A  summary  of  those  results  is  given  In  Appendix  A. 

The  working  equations  of  category  A  are  obtained  by  substituting 
Eqs.  (1)  and  (3)  for  pressure  into  the  equations  of  Appendix  A,  using 

o,  0,  k^  ■  constant,  8  ■  8(T),  k  •  <(T) 

wherever  these  quantities  occur  in  Eqs.  (1)  and  (3). 

Reduction  of  the  equations  to  forms  suitable  for  computation  fol¬ 
lows  a  pattern  very  close  to  that  of  the  test  problem  in  Ref.  1.  Those 
who  are  interested  in  pursuing  the  details  should  use  the  test  problem 
as  a  guide.  We  limit  ourselves  here  to  a  presentation  of  the  final 
results,  which  appear  in  Appendix  B. 

The  equations  of  category  B  follow  a  pattern  of  development  en¬ 
tirely  similar  to  those  of  category  A.  For  convenience,  we  have  taken 
k  to  be  a  function  of  x,  which  is  defined  by 
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-  .  i-1 
x  =  - 


l  nikl 


V(T  +  0)‘ 


(5) 


Requirements  for  the  equivalence  of  x  and  x  as  choices  of  indepen¬ 
dent  variable  for  x  follow  from  their  relation  to  one  another. 


xx  •  x  .  (6) 

The  values  of  the  variables  k ,  x,  and  x  are  all  greater  than  zero. 
Therefore,  for  each  continuous  functional  dependence  H(x)  assigned  to 
x,  there  is  a  reciprocal  one-to-one  correspondence  between  x  and  x  via 
the  relation  x/H(x)  *  x,  provided  the  function  x/H(x)  is  monotone  for 
positive  values  of  x  and  H(x) .  The  functions  we  will  use, 

x(x)  ■  Kx*,  K,y  •  constants  >  0, 

satisfy  these  requirements  on  x(x). 

The  equations  which  are  used  for  computations  based  on  parameters 
in  category  B  appear  in  Appendix  C. 


17 


V.  HUGONIOT  COMPUTATIONS.  REGION  I 


ALL  PARAMETERS  TREATED  AS  CONSTANTS 

Our  parametric  studies  are  performed  on  the  p-v  shock  Hugonlot 

-6  3  3 

centered  at  p  ■  10  megabars,  v  -  1  cm  /gm  (V  ■Hv  cm  ),  which 
o  o  000 

corresponds  to  ambient  fresh  water.  There  Is  no  loss  In  the  accuracy 
of  our  results  If  we  set  pQ  equal  to  zero.  The  Hugonlot  relation 
becomes 


E(T,  V,  a)  -  \  p(T,  V,  a)(Vo  -  V)  .  (7) 

The  Internal  energy  function  E(T,  V,  n)  has  been  defined  so  that  Eq  - 
0.  This  explains  the  absence  of  Eq  In  Eq.  (7). 

Our  first  Hugonlot  parameter  study  covers  Region  I,  and  Is  based 
entirely  on  constant  parameters.  For  this  particular  case  the  equa¬ 
tions  of  Appendix  B  and  those  of  Appendix  C  are  equivalent.  The  equa¬ 
tions  of  Appendix  B  are  moat  easily  reduced  to  the  present  case.  The 
computer  runs  made  with  constant  values  of  the  parameters  are  used  to 
decide  which  of  the  four  sets  of  covolume  factors,  k^,  appearing  in 
Table  4  should  be  used. 

We  use  Mader's^^  recommended  values  for  k,  B,  a,  and  0  based  on 
detonation  calculations  using  RDX  as  the  explosive.  They  have  the 
values  shown  in  set  1,  Table  5. 

We  chose  the  parameters  in  Ref.  14  based  on  the  explosive  RDX 
rather  than  TNT  because  our  system  contains  no  carbon  (see  Ref.  14).* 
The  method  of  computations  used  Is  discussed  briefly  in  Appendix  D. 
Computations  have  been  carried  out  in  Region  I  for  all  four  sets  of 
covolume  factors  to  determine  which  set  comes  closest  to  reproducing 
the  Thomas-Fermi  Hugonlot  at  high  pressures.  The  results  are  shown 
in  Fig.  3,  along  with  the  Hugonlot  results  of  Ref.  12,  using  an  older 

*The  four  sets  of  covolume  factors  in  Table  4  are  those  of  Ref. 
12,  with  covolime  factors  for  Hj,  Oj,  and  H^O  changed  to  the  values 
recommended  In  Ref.  14. 
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Tabic  4 

COVOLUME  FACTORS  AND  WEIGHTED  SUMS  OP  COVOLUME A FACTORS 
USED  IN  BKW  PARAMETER  STUDIES;  REGION  I 


Set  1 _ Set  2  Set  3  Set  4 


c 

c 

c 

c 

Species 

Number 

ki 

2  Vj 

J-i 

ki 

2vi 

j-i 

ki 

j-1 

ki 

IVj 

J-1 

1 

80 

-  20 

80 

-  20 

120 

-  30 

120 

-  30 

2 

65 

-  35 

65 

-  35 

98 

-  52 

98 

-  52 

3 

54 

-  46 

54 

-  46 

80 

-  70 

80 

-  70 

4 

45 

-  55 

45 

-  55 

67 

-  83 

67 

-  83 

5 

38 

-  62 

38 

-  62 

57 

-  93 

57 

-  93 

6 

0.3 

-100 

0.3 

-100 

0.5 

-150 

0.5 

-150 

7 

0.3 

-100 

0.3 

-100 

0.4 

-150 

0.4 

-150 

8 

0 

-100 

0 

-100 

0 

-150 

0 

-150 

9 

0 

-  20 

0 

-  40 

0 

-  30 

0 

-  60 

10 

350 

150 

350 

150 

350 

50 

350 

50 

11 

255 

55 

255 

55 

255 

-  45 

255 

-  45 

12 

180 

140 

180 

100 

180 

120 

180 

60 

13 

50 

10 

100 

20 

50 

-  10 

100 

-  20 

14 

250 

110 

250 

70 

250 

40 

250 

-  20 

15 

200 

80 

200 

60 

200 

20 

200 

-  10 

16 

230 

110 

230 

90 

230 

50 

230 

20 

17 

170 

50 

170 

30 

170 

-  10 

170 

-  40 

18 

115 

15 

115 

15 

173 

23 

173 

23 

19 

25 

5 

50 

10 

38 

8 

75 

15 

20 

240 

40 

240 

40 

240 

-  60 

240 

-  60 

21 

20 

40 

30 

60 

22 

100 

100 

150 

150 

23 

0 

0 

0 

0 

* 

All  covolumes  except  those  for  species  10,  12,  and  14  (H-0,  0,, 
H.O  respectively)  are  from  Ref.  12.  Covolumes  for  species  10,  12,  and 
14  are  from  Ref.  14. 


Table 
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Me  of  BKW  parameter*.  From  Pig.  3  we  that  all  four  solutions  based 
on  the  parameters  of  Raf.  14  and  the  (modified)  covolume  factors  of  Ref. 
12  represent  a  considerable  Improvement  at  the  lower  pressures  over  the 
solutions  appearing  in  Ref.  12.  However,  none  of  the  four  p-v  Hugonlota 
so  obtained  match  the  Thomas-Feral  Hugoniot  solutions  at  the  higher  pres¬ 
sures.  The  leftmost  Hugoniot  curve  (curve  1),  based  on  the  first  column 
of  covolume  factors  appearing  in  Table  4,  is  closest  to  the  Thomas-Feral 
Hugoniot.  We  therefore  chose  this  set  of  covoliae  factors  for  our  sub¬ 
sequent  programs. 


CONSTANT  k |  VARIABLE  S 

It  is  clear  from  Fig.  3  that  our  BXW  Hugoniot,  curve  (1),  is  too 
steep  at  the  high  pressures.  That  is,  for  the  portion  of  the  curve  in 
the  Thomas-Feral  region,  the  Hugoniot  pressures  corresponding  to  a  given 
specific  volume  are  much  higher  than  those  predicted  by  the  Thomas-Ferai 
model.  The  material  based  on  the  BKW  equation  with  constant  parameters 
appears,  therefore,  to  be  too  hard.  Our  first  attempts  tc  soften  the 
material  were  based  on  variations  of  S  with  temperature.  If  we  attribute 
the  steepness  of  the  Hugoniot  curve  to  excessive  departure  of  the  BKW 
equation  from  Ideality,  then  it  is  evident  from  Eqs.  (1)  and  (2)  and  the 
condition  x  2  0  that  we  want  8  to  decrease  as  the  Hugoniot  pressures  in¬ 
crease.  Since  temperatures  increase  on  the  Hugoniot  with  increasing 
pressures,  we  must  construct  functions  S(T)  that  have  the  appropriate 
value  of  0.16  at  a  temperature  of  5000*K,  and  that  decrease  monotonically 
with  temperature  for  temperatures  greater  than  5000*K. 

Two  functional  forms  meeting  these  conditions,  which  we  used  in 
Hugoniot  calculations,  are 


A  A 

8(T)  -  B2  + 


T 

5000 


40 


8 


2 


(8) 


constant , 


5000  STS  200 ,000*K  ; 
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(9) 

«  •  constant,  5000  S  T  £  10^*K  , 

A 

where  8(T)  =  8(T)/.16  and  T  is  measured  in  degrees  Kelvin. 

Equation  (8)  is  a  linear  temperature  relation  for  8  between  the 

A 

value  of  one  at  5000*K  to  an  arbitrary  value  8.  at  a  higher  temperature , 

*  A 

arbitrarily  chosen  as  200,000*K.  Hugoniot  runs  were  made  with  8^ 
taken  as  low  as  0.1.  The  asymptotic  approach  of  8  to  zero  in  Eq.  (9), 

T  •*  «*,  8  -*■  0,  can  be  hastened  or  retarded  by  appropriate  choices  of  the 
constant  a.  We  have  made  Hugoniot  calculations  using  Eq.  (9)  with  val¬ 
ues  of  a  between  0.01  and  4.5.  None  of  the  p-v  Hugoniot  curves  where 
Eqs.  (8)  and  (9)  were  used  exhibited  sufficient  departures  from. the 

A 

corresponding  Hugoniot  based  on  6  *1  to  account  for  the  discrepancy 
between  Thomas-Ferml  and  BKW  Hugoniot  results.  It  thus  appeared  that 
suitable  alterations  of  curve  (1),  Fig.  3,  could  not  be  made  by  varying 
the  parameter  8  over  a  reasonable  range  of  values  in  Region  II.  This 
observation  has  led  us  to  set  8  equal  to  0.16  in  the  remainder  of  our 
Investigations,  and  to  focus  our  attention  on  constructing  suitable 
functional  forms  for  the  parameter  x  to  achieve  a  smooth  connection 
between  the  experimental  and  Thomas-Ferml  Hugonlots  of  Fig.  2. 

x  AS  A  FUNCTION  OF  TEMPERATURE 

Another  possibility  for  obtaining  smaller  fluid  volumes  at  a  given 
pressure  than  those  given  by  our  constant-parameter  Hugoniot  is  to  allow 
the  covolume  factors  to  decrease  with  increasing  pressure  on  the  Hugo¬ 
niot.  Although  this  effect  does  occur  already  in  our  model  via  the  tem¬ 
perature  dependence  of  x,  the  covolume  factors  can  be  further  decreased 
by  varying  x,  which  appears  as  one  of  the  products  of  the  covolume  fac¬ 
tor  terms  in  the  definition  of  x. 

Two  possible  variations  of  x  have  already  been  discussed,  i.e., 
k(T)  and  x(x).  We  will  first  consider  variations  of  x  with  temperature. 

In  order  to  make  the  material  softer,  x  must  decrease  at  the  high¬ 
er  pressures.  Since  temperature  increases  with  Increasing  Hugoniot 
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pressures,  we  oust  construct  functions  of  temperature  for  k,  which  give 
the  correct  value  of  10.91  at  a  temperature  of  5000*K,  and  decreaae  (but 
remain  positive)  for  temperatures  greater  than  5000*K.  An  additional 
desirable  condition  would  be  to  have  the  temperature  derivative,  tc'(T), 
vanish  at  T  •  5000*K.  The  desirability  of  this  follows  from  the  appear¬ 
ance  of  k'(T)  in  the  equations  of  Appendix  B.  This  condition  has  been 
relaxed  because  it  had  no  significant  effect  on  Hugonlot  calculations 
on  or  near  5000*K.  Computations  were  made  with  two  different  functional 
forms  for  k(T)  ,  both  of  which  had  nonzero  values  for  k'(5000*K).  The  p-v 
Hugonlot  solutions  for  these  cases  were  essentially  identical  with  the 
p-v  Hugonlot  based  on  constant  k,  for  temperatures  between  5000*K  and 
15,000*K. 

Two  functional  forms  of  k(T)  have  been  chosen.  The  normalized 
forms  of  these  equations  are 


k<T)  -  <2  ♦  ^  5§00  “  > 


*2  "  constant,  5000  s  T  i  200,000*K  ; 


(10) 


tc(T)  -  1  + 


/T  -  5000  N"" 

V  5000  A  /  » 


A,  m  •  constants,  5000  S  T  S  10  *K  , 


6., 


(ID 


where  <(T)  s  <(T)/10.91,  and  T  is  again  measured  In  degrees  Kelvin. 

Equation  (10)  has  the  same  form  as  Eq.  (8),  and  the  parameter 
may  be  interpreted  in  the  same  way;  It  is  the  value  of  k  at  a  tempera¬ 
ture  of  200,000*K.  The  primary  reason  for  using  Eq.  (10)  for  ie(T)  was 
to  ascertain  whether  reasonable  variations  in  k  could  be  made  to  yield 
Hugonlot  curves  roughly  in  agreement  with  the  Thomas-Feral  Hugonlot  at 
pressures  of  approximately  20  megabars.  Results  of  calculations  using 
Eq.  (10)  are  shown  in  Fig.  4.  These  results  show  that  It  is  possible 
to  obtain  fair  agreement  with  the  Thomas-Feral  Hugonlot  at  the  lower 
end  of  the  Thomas-Fermi  region,  and  that  the  value  of  the  parameter  ie 


p  (megaban) 


Fig. 4—  BKW  Hugoniot 
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at  200,000*K  muse  be  reduced  Co  approximately  40  percent  of  the  value 

A 

at  5000*K.  The  linear  relation  for  <(T)  with  <^  3  0.4  predicts  negative 
values  for  ic  at  temperatures  greater  than  approximately  330,000*K.  This 
undesirable  feature  of  Eq.  (10)  is  not  shared  by  Eq.  (11),  which  was 
used  for  further  computations. 

The  curves  appearing  in  Fig.  5  were  computed  to  obtain  an  estimate 
of  the  correct  values  that  should  be  used  for  the  parameters  m  and  A  in 
Eq.  (11).  The  straight  line  in  this  figure  represents  the  linear  rela- 

A  A 

Cion ,  Eq.  (10),  for  <(T),  with  a  value  for  of  .40.  This  particular 
value  was  chosen  because  it  gave  the  best  agreement  between  the  Thomas- 
Fermi  and  BKW  Hugonlots  in  Fig.  4.  The  remaining  curves  in  Fig.  5  pro¬ 
vide  a  sample  of  <(T)  relations  obtained  from  Eq.  (11),  using  various 
values  f„r  the  parameters  m  and  A.  Approximate  limits  on  the  variation 
of  A  for  a  given  m  were  obtained  by  requiring  all  eligible  curves  to  go 

A 

through  a  band,  .40(1  -  S^)  2  <(200,000)  2  .40(1  +  S^) ,  where  the  values 
of  and  are  arbitrarily  1/2  and  1,  respectively.  Using  this  method 
and  favoring  those  curves  which  lay  nearer  to  the  linear  relation  for 
temperatures  between  5000*K  and  200,000*K,  we  have  chosen  a  value  of  1 
for  m,  and  a  set  of  A's  shown  in  Fig.  6. 

Our  final  choice  of  parameters  for  <(T)  is  A  ■  15,  m  ■  1,  based  on 
results  shown  in  Fig.  6.  The  functional  form  of  <(T)  which  we  therefore 
accepted  as  final  is 


*(T) 


10.91 
T  -  5000 
75,000 


(12) 


All  remaining  calculations,  based  on  <(T)  both  on  and  off  the  Hugoniot 
in  Region  I,  will  use  this  equation.  Calculations  off  the  Hugoniot 
are  discussed  in  a  later  section. 


<  AS  A  FUNCTION  OF  x 

Before  we  can  construct  functions  <(x)  that  decrease  with  increas¬ 
ing  Hugoniot  pressures  and  temperatures,  we  must  have  some  knowledge 
of  how  x  varies  on  the  Hugoniot.  No  clear  statement  of  x's  variation 
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on  the  Hugonlot  can  be  made  directly  from  its  definition, 

23 

x  -  ~ - -  , 

V(T  +  8)° 

for  while  V  decreases  with  increasing  p,  (T  +  6)°  increases,  and  the 
variations  of  the  nunerator  are  unpredictable  without  extensive  calcu¬ 
lation. 

Some  insight  into  x's  variation  can  be  obtained  from  the  results 
of  Hugonlot  calculations  shown  in  Table  5.  These  results  are  based  on 
constant  values  for  all  of  the  parameters  (Mader's  recommended  values 
for  a,  B,  k,  6,  and  the  first  set  of  covolume  factors  in  Table  4).  In 
Table  5,  x  on  the  Hugonlot  decreases  monotonically  toward  zero  with  in¬ 
creasing  pressure  (and  temperature) .  Because  these  results  are  based 
on  constant  k,  they  are  only  an  indication  of  the  behavior  of  x  for 
variable  k.  However,  they  are  sufficient  to  establish  reasonable  trial 
forms  for  k(x).  Trial  forms  are  therefore  chosen  to  vanish  at  x  »  0, 
to  give  the  correct  value  of  k  -  10.91  at  T  »  5000°K  on  the  Hugonlot 
(where  x  «  0.50566  from  Table  6),  and  to  increase  monotonically  with 
x  for  x  >  0. 

Of  the  various  forms  tried,  the  following  simple  expression  most 
successfully  matches  Hugonlot  curves  to  the  Thomas-Fermi  Hugonlot  at 
high  pressures: 


"  (o. 50566)  ’  <13) 

x  >  0,  y  ■  constant  , 

where  k(x)  =  x(x)/10.91. 

Calculations  based  on  Eq.  (13)  used  in  association  with  the  rela¬ 
tions  of  Appendix  C  are  presented  in  Fig.  7,  which  is  a  set  of  p-v 
Hugoniot  curves  parameterized  on  y.  From  Fig.  7,  we  have  chosen  the 
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Table  6 


HUGONIOT  CALCULATIONS  BASED  ON  FIRST  SET  OF  PARAMETERS,  TABLE  5 
COVOLUME  FACTORS  FROM  SET  1  OF  TABLE  4 


T(°K)  x  10'3 

p(megabars) 

X 

5 

0.887 

0.50566 

7 

1.23 

0.475 

10 

1.60 

0.425 

15 

2.23 

0.365 

20 

2.85 

0.311 

30 

4.50 

0.205 

50 

6.50 

0.139 

70 

8.39 

0.113 

100 

11.34 

0.0909 

150 

17.15 

0.0698 

200 

24.12 

0.0574 

300 

40.15 

0.0420 

500 

78.20 

0.0270 

700 

119.1 

0.0196 

1000 

189.4 

0.0132 
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VI.  RESULTS  ON  ISOTHERMS  AND  ISOMETRICS.  REGION  I 

Table  5  contains  a  summary  of  the  parameters  used  for  the  three 
sets  of  computations  [k  -  constant.  <  ■  k(T),  k  ■  k(x)]  in  the  remainder 
of  this  Report. 

Pressure-volume-temperature  calculations  for  a  set  of  isotherms 
spanning  the  low-pressure  part  of  the  Thomas-Fermi  region  are  compared 
to  the  corresponding  Thomas-Fermi  isotherms  in  Figs.  8.9.  and  10,  where 
each  figure  employs  a  different  set  of  the  parameters  appearing  in  Table 
5.  Energy-volume-temperature  comparisons  are  made  in  Figs.  11,  12,  and 
13  on  the  same  set  of  Isotherms  that  appear  in  the  p-v-T  diagram. 

The  p-v-T  comparisons  of  Figs.  8,  9,  and  10  are  somewhat  surprising, 
because  the  solutions  for  the  constant-K  set  of  parameters,  which  gave 
the  worst  p-v  Hugoniot  results,  match  the  p-v  isotherms  of  the  Thomas- 
Fermi  model  better  than  do  the  variable-K  cases,  whose  pressures  on  iso¬ 
therms  are  consistently  lower  than  the  corresponding  Thomas-Fermi  pres¬ 
sures.  However,  the  e-v-T  comparisons  shown  in  Figs.  11,  12,  and  13 
yield  consistently  low  BKW  energies  on  Isotherms  for  all  three  sets  of 
parameters.  Comparison  of  p-v-T  and  e-v-T  results  for  each  set  of  cal¬ 
culations  suggests  that  the  variable-K  cases  might  yield  better  p-v-e 
results  than  the  constant-K  case,  since  such  results  are  the  temperature 
ellmlnant  of  the  p-v-T  and  e-v-T  calculations  in  Figs.  8  to  10  and  11 
to  13,  respectively.  This  does  indeed  turn  out  to  be  true,  as  can  be 
seen  in  Figs.  14,  15,  and  16.  In  these  figures,  the  BKW  p-e  curves 
corresponding  to  four  isometrics  spanning  the  low-pressure  Thomas-Fermi 
region  come  much  closer  to  the  corresponding  Thomas-Fermi  isometrics 
for  the  varlable-K  calculations.  This  last  set  of  results  is  consistent 
with  the  p-v  Hugoniot  results  presented  earlier,  for  p-v  Hugoniot  cal¬ 
culations  use  only  the  pressure-volume-energy  form  of  the  equation  of 
state.  Thus,  the  variable  k's  chosen  to  fit  known  p-v  Hugoniot  results 
were  tailored  to  favor  a  correct  form  of  the  p-v-e  relation. 

Since  our  equation  of  state  will  be  primarily  used  in  hydrodynamic 
calculations,  which  mainly  use  the  p-v-e  form,  we  base  our  final  choice 
of  parameters  on  the  comparisons  appearing  in  Figs.  14  to  16.  In  these 
figures,  the  k(T)  and  k(x)  calculations  yield  essentially  the  same 


p  (megabars) 


v  (cc/gm) 


Fig.  10— Thomas-Fermi  and  BKW  p-v  isotherms.  Region  I.  BKW  results 
are  based  on«(x)  given  by  Eq.  (14). 


Fig.  14— Thomas-Fermi  and  BKW  p-e  isometrics.  Region  I. 
BKW  results  are  based  on  tt=  10.91. 


p  (megabars) 
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Fig.  15— Thomas-Fermi  and  BKW  p-e  isometrics.  Region  I.  BKW  results 
are  based  on  k  (T)  given  by  Eq.  (12). 


p  (megabars) 
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result.  Both  predict  pressures  that  are  consistently  higher,  by  approx¬ 
imately  7  percent,  than  the  Thomas-Ferml  solution  at  the  same  energy 
and  volume. 

We  choose  the  k(T)  parameters  (set  2  of  Table  5)  to  represent  the 
p-v-e  results  of  this  Report  in  Region  I,  since  the  k(T)  calculations 
on  isotherms  agreed  a  little  better  than  the  k(x)  calculations  with 
Thomas-Ferml  results. 

Because  of  the  generally  poor  performance  of  the  variable-K  iso¬ 
thermal  calculations  and  the  equally  poor  constant-K  p-v-e  (and  Hugo- 
niot)  results,  we  will  not  present  detailed  p-v-T  or  even  e-v-T  results 
in  Regions  I  and  II  based  on  the  BXW  equation  or  any  of  its  modifica¬ 
tions  considered  here.  Detailed  pressure-volume-energy  results  are 
presented  later  in  Fig.  21. 
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VII.  HUGONIOT  CALCULATIONS.  REGION  II 

Hugonlot  calculations  with  all  three  secs  of  parameters  appearing 
in  Table  6  showed  chat  at  5000*K  and  7000*K  there  were  insignificant 
traces  of  all  species  other  than  HjO.  Thus,  below  5000*K,  the  Hugonlot 
calculations  were  based  on  pure  water.  This  is  not  true  however,  for 
p-v  points  off  of  the  Hugonlot.  Off  of  the  Hugonlot,  the  six  species 
H,  0,  H^,  0 HjO,  and  OH  were  assumed  to  exist  in  chemical  equilibrium 
in  Region  II.  Another  simplification  occurring  in  Region  II  stems  from 
the  fact  that  its  pressure  range,  which  runs  from  about  0.5  to  1  mega¬ 
bar,  approximately  coincides  with  a  large  part  of  the  range  of  pressures 
over  which  the  BKW  parameters  of  Ref.  14  were  established.  Hugonlot 
calculations  were  made  using  these  parameters  (set  1,  Table  5,  and  the 

k  of  set  1,  Table  4) ,  in  association  with  standard  state  data  given  at 
*  * 

100*K  increments  from  100*K  to  60Q0*K.  The  p-v  Hugonlot  curve  obtained 

from  these  calculations  is  placed  on  Fig.  17,  which  is  a  partial  repro¬ 
duction  of  Fig.  1,  Ref.  14.  The  results  of  our  Hugonlot  calculations 
can  be  seen  to  be  in  good  agreement  with  those  of  Ref.  14  above  approx¬ 
imately  .2  megabar,  with  minor  discrepancies  that  can  be  attributed  to 
differences  in  standard  state  data  and  ambient  temperatures  at  the  foot 
of  the  Hugonlot  between  our  calculations  and  those  of  Ref.  14.  Below 
approximately  .17  megabar  the  BKW  equation  predicts  shock  temperatures 
lower  than  the  ambient  temperature  of  the  unshocked  material.  This 
renders  comparisons  between  our  results  and  those  of  Ref.  14  somewhat 
meaningless  below  .17  megabar,  since  in  this  region  the  results  of 
either  set  of  calculations  represent  physical  nonsense. 

Although  the  p-v  Hugonlot  results  using  the  BKW  method  exhibit 
fairly  good  agreement  with  experimental  data,  the  p-T  Hugonlot  calcula¬ 
tions  lead  to  results  that  in  our  opinion  do  not. bear  even  a  casual  re¬ 
lation  to  the  p-T  Hugonlot  data  for  water.  This  observation  is  based 
on  a  comparison  of  the  experimental  and  BKW  p-T  Hugonlot  curves  shown 
in  Fig.  18.  The  results  of  Fig.  18  place  considerable  doubt  on  the 

+See  Ref.  24  for  standard  state  data  below  6000*K. 


p  (megabars) 


■44- 


Fig.  17— BKW  Hugoniot  curves.  Regions  il  and  III.  Solid  and  dashed  curves 
are  based  on  first  set  of  parameters  (Table  6),  with  =  250. 


Fig.  18— -Comparison  of  BKW  p-THugoniot  curve.  Regions  II  ond  III 
with  that  of  experimental  data. 
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value  of  Che  temperature-dependent  BKW  equations  for  representing  the 
equilibrium  properties  of  water  near  the  Hugonlot  at  pressures  above 
.1  megabar.  This  conclusion  Is  supported  by  Che  results  of  the  follow¬ 
ing  section. 
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VIII.  RESULTS  ON  ISOTHERMS  AND  ISOMETRICS.  REGION  II 


In  Region  II  it  is  difficult  to  obtain  a  good  evaluation  of  BKW 
results  on  isotherms  and  isometrics  at  pressures  above  .25  megabar  and 
for  temperatures  higher  than  a  few  thousand  degrees  Kelvin,  since  ref¬ 
erence  curves  are  difficult  to  obtain  in  these  regions.  The  best  that 

we  can  do  is  to  extrapolate  the  results  of  Rice  and  Walsh  to  higher 

(4) 

temperatures  and  pressures  than  were  incorporated  in  their  work. 

For  the  p-v  isotherms,  we  use  the  results  of  Ref.  8  as  reference  curves 
These  results,  however,  represent  an  extrapolation  of  one  of  the  assump 
tions  of  Ref.  4,  (Sv/ST)^  ■  f(p),  to  temperatures  as  high  as  10,000*K. 
No  justification  of  this  assumption  appears  in  the  text  of  Ref.  8.  For 
p-e  comparisons  on  isometrics,  we  obtain  reference  curves  to  .450  mega¬ 
bar  by  extrapolating  the  curve  fit  of  Ref.  4  for  the  quantity  (3h/3v) 

P 

£(p),  to  pressures  of  .450  megabar.  In  Ref.  4,  this  expression  is  ter¬ 
minated  at  a  maximum  pressure  of  .250  megabar.  The  rather  uncertain 
value  of  these  extrapolations  leads  us  to  make  comparisons  only  in  a 
small  neighborhood  of  the  Hugoniot  curve  above  .250  megabar,  since  the 
extrapolations  are  more  reliable  near  the  Hugoniot  than  anywhere  else. 

In  Region  II  as  in  Region  I,  the  BKW  equation  exhibits  a  poor  tem¬ 
perature  performance,  as  is  clearly  borne  out  by  Fig.  19.  This  figure 
compares  the  BKW  p-v-T  solutions  based  on  the  first  set  of  parameters 
of  Table  5  with  the  results  of  Ref.  8.  These  results,  along  with  equal 
ly  unflattering  comparisons  of  the  BKW  e-v  isotherms  with  reference 
curves  (not  shown)  has  led  us  to  conclude  that  in  Region  II  as  well  as 
in  Region  I,  the  BKW  method  does  not  yield  acceptable  p-v-T  and  e-v-T 

4. 

equations  of  state  near  the  Hugoniot  curve.’ 

Pressure-energy  isometrics  are  shown  in  Fig.  20  for  both  the  ref¬ 
erence  and  BKW  systems.  A  comparison  of  the  reference  and  BKW  isomet¬ 
rics  of  Fig.  20  has  led  us  to  reject  the  BKW  p-v-e  results  at  pressures 
lower  than  approximately  .3  megabar. 


Reference  isotherms  and  Isometrics  havu  also  been  obtained  from 
the  more  recent  results  of  Ref.  25.  Comparisons  of  the  BKW  curves  of 
Figs.  19  and  20  with  these  later  results  have  not  produced  any  changes 
in  the  conclusions  and  recommendations  of  this  Report. 
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IX.  SUMMARY  OF  RESULTS  AND  DISCUSSION 


1.  Our  studies  have  shown  that  the  unmodified  forms  of  the  BKW 
equation  are  generally  unsatisfactory  near  the  Hugoniot.  In  particular, 
the  p-T  Hugoniot  curve  predicted  by  the  BKW  method  is  completely  in 
error.  However,  results  on  p-v  isotherms  above  200, 000 *K  suggest  that 
the  unmodified  form  of  the  BKW  p-v-T  relation  is  satisfactory  above 
this  temperature  if  the  parameters  of  Ref.  13  are  used. 

2.  Of  the  modifications  attempted,  only  those  that  altered  the 
covolune  factors  of  the  mixture  showed  any  appreciable  tendency  to 
bring  the  p-v  Hugoniot  back  into  line  with  the  known  Thomas-Fermi  re¬ 
sults.  Both  of  the  trial  forms  for  variable  covolume  factor  [k(x)  and 
k(T)]  produced  satisfactory  p-v  Hugoniot  results.  In  each  case  this 
was  accomplished  by  forcing  the  modified  forms  of  the  BKW  equation  to 
lie  closer  to  the  ideal-gas  equation  of  state  than  does  the  unmodified 
form. 

3.  The  modified  BKW  forms  produced  unsatisfactory  p-v-T  and  e-v-T 
results  at  temperatures  as  high  as  7Q0,000*K,  except  in  the  ideal-gas 
range.  However,  the  modified  forms  (in  particular  the  form  where  k 
assumed  a  dependence  upon  temperature)  produced  reasonable  p-v-e  re¬ 
sults  above  approximately  300  kilobars  for  the  region  lying  between 
the  Hugoniot  curve  and  the  10,000  megabar-adiabat . 

4.  The  combined  p-v-e  results  of  Regions  I  and  II  are  presented 
in  Fig.  21.  Results  have  been  altered  slightly  so  that  they  smoothly 
Join  the  Thomas-Fermi  solutions  of  Ref.  2,  which  also  appear  in  Fig. 

21.  The  dashed  lines  represent  linear  extrapolations  of  computed  iso¬ 
metrics  (solid  lines) .  The  figure  also  includes  the  Hugoniot  curve 
and  the  10,000-megabar  adiabat  (the  isentrope  through  the  Hugoniot  at 

4 

a  pressure  of  10  megabars) .  These  curves  represent  the  approximate 
boundaries  of  the  region  required  for  nuclear  burst  calculations  in 
water. ^  The  10,000-megabar  adiabat  was  obtained  by  numerically  inte¬ 
grating  the  relation  de  -  -pdv,  with  the  aid  of  the  p-v-e  results  of 
Fig.  2i,  from  the  10,000-megabar  point  on  the  Hugoniot. 


See  Ref,  1  for  a  discussion  of  the  approximate  range  requir  ments 
for  nuclear  burst  calculations. 


Hugoniot 


Fig. 21— Pressure-volume-energy  relation  for  water  obove  300  kilometers 
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W«  regard  the  p-v-a  results  obtained  hare  aa  somewhat  unsatisfactory, 
but  perhaps  as  good  aa  could  reasonably  be  provided  by  the  BXK  fora  of 
the  thermal  equation  of  state.  Any  further  improvement  in  the  combined 
p-v-T,  e-v-T,  and  p-v-e  relations  will  simply  require  a  more  satisfactory 
point  of  departure  (thermal  equation  of  state,  partition  function,  etc.) 
than  is  offered  by  the  BKW  equation. 
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(A-9) 


(A- 10) 


(A-ll) 


The  following  definiciona  apply  Co  Che  equations  of  this  summary: 

T  -  298. 16®K 
o 

o  .  , 

p  -  1  atm 


,,o  ,  o 

V  -  ngKT/p 


volume  of  n  mols  (M  go)  of  mixture  In  the 
standard  state  ° 


mQ  £  J>| 

H^(T)  ”  molar  enthalpy  of  the  1  species  in  the  standard  state 
at  temperature  T  relative  to  the  enthalpy  of  the  refer¬ 
ence  materials  Hj,  0^,  and  e  in  the  standard  state  at 
temperature 

H°(Tq)  ■  enthalpy  of  formation  of  one  mol  of  the  ith  species  in 
its  standard  state  from  the  reference  materials  (0_, 

H2,  e)  in  their  standard  states  at  temperature  Tq 

F?(T)  •  molar  Gibbs  free  energy  of  the  iC^  species  in  its 
standard  state  at  temperature  T 


in  K  (T) 
pi 


V 


ij 


c  v. .F?(T) 


-  I  ^ 

j-1  RT 


stoichiometric  coefficient  of  the  component  ap¬ 
pearing  in  the  i  reaction 
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y ^  -  sola  of  atomic  hydrogen  in  Mq  gm  of  mixture 

?2  “  n22  “  ■ol*  of  atomic  oxygen  in  Mq  gm  of  mixture 

^3  "  n23  "  ao^*  electrons  in  Mq  gm  of  mix turn 

The  standard  state  data  for  fifteen  temperatures  from  5G00*K  to 
1,000,000*K  (Region  1)  are  tabulated  in  Ref.  11.  Host  of  these  data 
can  also  be  found  in  Ref.  12.  The  standard  Btate  data  in  this  region 
are  given  at  the  following  temperatures  (in  degrees  Kelvin):  5000, 

7000,  10,000,  15,000,  20,000,  30,000,  50,000,  70,000,  100,000,  150,000, 
200,000,  300,000,  500,000,  700,000,  1,000,000. 

These  data  have  been  used  for  all  diatomic  species  and  for  U^O. 
However,  all  standard  state  data  for  tne  monatomic  species  above 
5000*K  have  been  obtained  from  a  partition-function  program  developed 
by  P.  Gilmore  at  Rand.  These  (unpublished)  results  have  been  incor¬ 
porated  because  they  represent  more  recent  data,  and  yield  a  signifi¬ 
cantly  closer  matching  with  the  standard-state  data  of  Ref.  24,  which 
is  used  for  the  lower  temperatures. 

Standard  state  data  for  Region  II  are  taken  from  the  JANAP  ther- 
(24) 

mochemical  tables.  The  data  are  given  for  temperatures  from  0 

to  6000*K  in  100-degree  increments,  plus  data  at  298*K. 

Stoichiometric  coefficients  and  r^  appear  in  Tables  2  and 
3,  respectively. 
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Appendlx  E 

SUMMARY  OP  THERMODYNAMIC  RELATIONS  BASED  ON  THE  MODIFIED 
BKW  EQUATION  WITH  0  -  l(T),  H  -  k(T) 

The  working  equations  for  Regions  I  and  II,  based  on  the  modified 
BKW  equation  with  variable  parameters,  0  -  0(T),  k  -  k(T),  are  pre¬ 
sented  here.  These  results  were  obtained  by  placing  into  the  equations 
of  Appendix  A  the  thermal  equation  of  state  [p  -  p(T,  V,  n)]  shown  be¬ 
low,  and  carrying  out  the  differentiation  and  integrations  wherever 
it  was  possible  to  do  so.  The  equations  are: 


*  - - — - 

V(T  +  9) 

23 

k  -  h(T)  £ 
i-1 

23 

n  *  )  n. 
g  Li 

i-1 


(B-l) 

(B-2) 

(B-3) 

(B-4) 


23  23 

E  .  I  n^d)  -  !?<*„)]  +  l  «,«Jd0> 


(B-5) 


i-1 


i-1 


+  68317  -  V1  +  Vt2{[it-0  - 


0Z(T)  L 


the  acthalpy  of  formation  of  electrona,  H^V'  **r°’  ***  th*lr 
molar  enthalpy  difference  ia  glvan  by  tha  relation 

-  »L<V  •  f  »(T '  V 


r,.  r„,  -r 


_  *  V  2^V  1  <  J  S  20 

J  Jyl  72  y3  * 


(B-6) 


o.^i-1 


i  . «  (i>  (£9 

J  pj  xaT y 


rsr  >l£!*L iil 

*  IQ  vi^ - 5(T) 


+  fY'(T)x  IvJ 


(B-7) 


'l+  1  ‘u”l* 


(B-8) 


’i+i  r2i°t  • 


(B-9) 


■  I  r3i”i  ' 


(B-10) 


v  “  V/M 


(B-ll) 


e  -  E/M 


(B-12) 


Sea  Appendix  A  for  definition*  of  the  symbol*  used  in  this  summery  of 

equations.  Values  for  the  stoichiometric  coefficient  v. .  and  r  ,  can 

i  j  on 

be  found  in  Tables  2  and  3.  Table  3  also  Includes  the  sum 

23 

Im 

J-i 

The  covolume  factors  and  the  sums 

23 

I  VJikj 
j-1 


appear  in  Table  4.  Sources  for  the  necessary  standard  state  data  are 
given  in  Appendix  A. 
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Appendlx  C 

SUMMARY  OP  THERMODYNAMIC  RELATIONS  BASED  ON  THE 
MODIFIED  BKW  EQUATION  WITH  B  -  B(T),  k  -  k(x) 

The  working  equations  for  Regions  I  and  II,  based  on  Che  modified 
BKW  equation  with  variable  parameters  B  -  B(T) ,  k  -  k(x) ,  are  pre¬ 
sented  here.  These  results  were  obtained  by  placing  into  the  equa¬ 
tions  of  Appendix  A  the  theraal  equation  of  state  [p  -  p(T,  V,  n)] 
shown  below,  and  carrying  out  the  differentiations  and  integrations 
wherever  it  was  possible  to  do  so.  The  equations  are: 


n  Kt  .... 

p  -  u  + 

x  -  xx(x) 

i — l— 

V(T  +  0)“ 

23 

£-  l  ki°i 

i-1 

23 

n  -  )  n, 

g  L  i 

i«l 


(C-l) 

(C-2) 

(C-3) 

(C-4) 

(C-5) 


The  enthalpy  of  formation  of  electrons,  H23(t0).  t*  *«<>.  •«* 
their  molar  enthalpy  difference  la  given  by  the  relation 


H^(I)  -  H^(I0)  -  f  R(t  -  T0> 


lJ  "  KJ,l1J,221y3  31 '  1  s  S  s 


^  ■  vcoogr1 11  -  -  [ c f  ■ o  t  «  * 

j  ° 


♦*(I  VO 

^  -f _ 1 


'l  +  I  rU°l  ■ 


y2  +  I  r2l“l 


-  1  r3tni‘ 


v  -  V/M 


(C-12 


.3 
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See  Appendix  A  for  definitions  of  the  symbols  used  in  this  stsomary  of 
equations.  Values  for  the  stoichiometric  coefficients  and  r^  can 
be  found  in  Tables  2  and  3.  Table  3  also  includes  the  sum 

23 

I  ’ij 

J-l 

The  covolume  factors  and  the  sums 

23 

Z  Wjiki 

J-l 

appear  in  Table  4.  Sources  for  the  necessary  standard  state  data  are 
given  in  Appendix  A. 


"  «  "W..  "VWV  .V-^SK,V>  . 
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Appendix  D 

BRIEF  OUTLINE  OF  THE  COMPUTATIONAL  METHOD  FOR 
EQUILIBRIUM  CALCULATIONS  IN  REGIONS  I  AND  II 

The  NewCon-Raphson  iCeraCion  method  is  applied  to  the  set  of 
five  simultaneous  nonlinear  equations  in  the  five  variables,  y^,  jy 
jy  y^,  jy  The  equations  are  solved  for  these  variables  as  functions 
of  volume  and  temperature.  For  the  Hugoniot  calculations,  a  sixth  re¬ 
lation  is  added  to  the  set  of  equations  and  the  volume  y^  is  added 
to  the  list  of  variables. 

The  five  equations  are  f^  -  0,  i  -  1  ...  5,  where  the  f^  are  de¬ 
fined  by  the  following  relations 

20 

h  *  rn“i  '  2 

i-1 

20 

£2  ‘  y2  +  1  r2l“i  '  1 
i-1 

20 

f3  ’  y3  -  I  r3i“i 
i-1 

20 

f4  “  y4  -  yl  '  y2  -  y3  *  I  °i 

i-1 

23 

f5  -  y5  -  H(T)  l  klBi 
i-1 

The  Hugoniot  relation  is  f^  -  0,  where 
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...  ,„iKf,r  .^p,.  >t| 


,  g  V0  I  y6 

6  3  2y6 

y48T(l  +  xe  x) 

Appendices  A,  B,  and  C  contain  definitions  of  symbols  used  in 
these  equations.  In  addition,  the  variables  y^,  y^,  and  y^  have  the 
following  definitions. 

y4  ■  ■  total  taols  of  mixture 

y^  ■  k  »  k(T)  [average  covoluae  factor  of  mixture] 

y^  ■  v  «  specific  volume 

Linear  estimates  of  the  above  equations,  which  are  used  in  the 
iterative  process,  can  be  represented  in  the  following  form: 

£i(yJ)+I(^-)  dyJ"° 

J  3 

where 


k+1 

J 


k  ,  k 
jj  +  lyj 


i,  J»l,  ...  ,5  for  non-Hugoniot  calculations 
1,  J«l,  ...  ,6  for  Hugoniot  calculations 

In  these  relations,  the  superscript  k  denotes  the  kC^  Iterate. 
For  calculations  on  the  Hugoniot,  the  Initial  guesses,  y^ ,  were  taken 
as  the  values  given  In  Ref.  12.  Initial  guesses  for  y^  along  an  Iso¬ 
therm  were  taken  from  results  of  calculation  of  neighboring  volume 
and  temperature  points. 
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Preaaure  and  energy  reaulca  can  be  obtained  directly  after  the 
determination  of  the  variablea  y^,  J«l,  ...  ,5  or  6,  at  a  given  T  and 
v  or  T,  reapectively. 

The  accuracy  criterion  for  convergence  of  the  iterative  proceaa 

la 


where  max  aignifiea  the  maximum  value  overall  j,  at  a  given  iteration 
(k).  This  criterion  alao  givea  the  reault 

lci<  io'8 

i 

Thla  program  ia  written  in  the  MAP  language  and  run  on  the 
IBM  7040/44  ayatem. 


-65- 


REFERENCES 


1.  Papetti,  R.  A.,  and  M.  Fujlsaki,  Thermoohemical  Relatione  for  Water 
at  High  Temperaturee  and  Preeeuree,  Tha  Rand  Corporation,  RM-4970-PR, 
Dacaabar  1966. 

2.  Lattar,  A.,  and  R.  Latter,  Equation  of  State  of  Water ,  The  Rand 
Corporation,  RM-1492-AEC,  May  1955. 

3.  Al'tschuler,  L.  V.,  A.  A.  Bakanova,  and  R.  F.  Trunin,  Phaae  Trans¬ 
formations  of  Vater  Compressed  by  Strong  Shock  Waves,  Soviet  Phyeice , 
Doklady,  Vol.  3,  No.  4,  January  1958,  pp.  761-763. 

4.  Rice,  M.  H.,  and  J.  M.  Walsh,  Equation  of  State  of  Water  to  250 
Kilobars,  J.  Chem.  Phye.,  Vol.  26,  No.  4,  April  1957,  pp.  824-830. 

5.  Feynman,  R.  P.,  N.  Metropolis,  and  E.  Teller,  Equations  of  State  of 
Elements  Based  on  the  Generalized  Fermi -Thomas  Theory,  Phye.  Rev. , 
Vol.  75,  No.  10,  May  1949,  pp.  1561-1572. 

6.  Knopoff,  L.,  Equation  of  State  for  Water  at  High  Pressure,  J.  Chem. 
Phye.,  Vol.  28,  No.  6,  J»ma  1958,  pp.  1067-1069. 

7.  Walker,  W.  A.,  and  H.  M.  Sternberg,  The  Chapman -Jouguet  Isentrope 
and  the  Underwater  Shockwave  Performance  of  Pentolite,  in  The  Fourth 
Syirpoeium  on  Detonation,  October  12-15,  1965,  Preprints,  Vol.  1, 

U.S.  Naval  Ordnance  Laboratory,  White  Oak,  Silver  Spring,  Maryland, 
pp.  156-169. 

8.  Howard,  J.  C.,  Thermodynamio  Data  for  Water,  University  of  Califor¬ 
nia  Radiation  Laboratory,  UCRL-6455-T,  April  1961. 

9.  Snay,  H.  G. ,  and  1.  Stegun,  Auxiliary  Funotione  for  Thermodynamic 
Calculations  Based  on  the  Kietiakoueky-Wilson  Equation  of  State 
(Rirkwood-Montroll  I  -  Funotione) ,  U.S.  Naval  Ordnance  Laboratory, 
NAVORD  1732,  January  1951. 

10.  Rosenbaum,  J.  H.,  An  Equation  of  State  for  Water  at  Extreme  Pree- 
euree,  U.S.  Naval  Ordnance  Laboratory,  NAVORD  3847,  November  1954. 

11.  Snay,  H.  G,,  and  J.  F.  Butler,  Equation  of  State  for  Water ,  U.S. 
Naval  Ordnance  Laboratory,  NAVORD  4181,  November  1956. 

12.  Gleysal,  A.  N.,  and  H.  G.  Snay,  Equation  of  State  of  Water,  U.S. 
Naval  Ordnance  Laboratory,  NAVORD  6749,  December  1959. 

13.  Mader,  C.  L.,  Detonation  Performance  Calculations  Using  the 
Kietiakou8ky-Wil8on  Equation  of  State,  Los  Alamos  Scientific  Lab¬ 
oratory,  LA  2613,  October  1961. 


-66- 


14.  Mader,  C.  L.,  Detonation  Properties  of  Condensed  Explosives  Com¬ 
puted  Using  the  Beoker-Kistiakowsky -Wilson  Equation  of  State,  Lot 
Alamos  Scientific  Laboratory,  LA  2900,  February  1963. 

15.  Covan,  R.  D.,  and  W.  Fickett,  Calculation  of  the  Detonation  Prop¬ 
erties  of  Solid  Explosives  with  the  Kistiakowsky-Wilaon  Equation 
of  State,  J.  Chem.  Phys.,  Vol.  24,  No.  5,  Hay  1956,  pp.  932-939. 

16.  Christian,  E.  A.,  and  H.  G.  Snay,  Analysis  of  Experimental  Data 
on  Detonation  Velocities ,  U.S.  Naval  Ordnance  Laboratory,  NAVORD 
1508,  February  1951. 

17.  Courant,  R. ,  and  K.  0.  Friedrichs,  Supersonic  Flou  and  Shook  Waves, 
Interscience  Publishers,  Inc.,  New  York,  1948. 

18.  Fickett,  W.,  Intermoleoular  Potential  Funotions  for  Some  Simple 
Molecules  From  Available  Experimental  Data,  Los  Alamos  Scientific 
Laboratory,  LA  2665,  July  1962. 

19.  Fickett,  W.,  Calculation  of  the  Detonation  Properties  of  Condensed 
Explosives,  Physios  of  Fluids ,  Vol.  6,  No.  7,  July  1963,  pp.  997- 
1006. 

20.  Fickett,  W.,  W.  W.  Wood,  and  Z.  W.  Salsburg,  Investigations  of  the 
Detonation  Properties  of  Condensed  Explosives  with  Equations  of 
State  Based  on  Intermolecular  Potentials.  I.  RDX  With  Fixed  Product 
Composition,  J.  Chem.  Phys.,  Vol.  27,  No.  6,  December  1957,  pp.  1324- 
1329. 

21.  Fickett,  W. ,  Detonation  Properties  of  Condensed  Explosives  Calcu¬ 
lated  with  an  Equation  of  State  Based  on  Intermolecular  Potentials, 
Los  Alamos  Scientific  Laboratory,  LA  2712,  May  1962. 

22.  Keeler,  R.  N.,  B.  J.  Alder,  and  M.  van  Thiel,  Corresponding  States 
at  Small  Interatomic  Distances,  University  of  California  Radiation 
Laboratory,  UCRL-7907  Rev.  I,  November  1964. 

23.  Hlrschfelder ,  J.  0.,  C.  F.  Curtiss,  and  R.  B.  Bird,  Molecular 
Theory  of  Gases  and  Liquids,  Wiley,  New  York,  1954. 

24.  Joint  Army-Navy-Air  Force  Thermochemical  Panel,  JANAF  Interim 
Thermoahemical  Tables,  (H.0) ,  prepared  by  the  Thermal  Laboratory 
of  the  Dow  Chemical  Company,  March  1961. 

25.  Papetti,  R.  A.,  and  M.  Fujisaki,  The  Rice  and  Walsh  Equation  of 
State  for  Water:  Discussion,  Limitations,  and  Extension,  J.  Appl. 
Phys.,  Vol.  39,  No.  12,  November  1968,  pp.  5412-5421. 


